Minimization of nonresonant effects in a scalable Ising spin 
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The errors caused by the transitions with large frequency offsets (nonresonant 
transitions) are calculated analytically for a scalable solid-state quantum computer 
based on a one-dimensional spin chain with Ising interactions between neighboring 



spins. Selective excitations of the spins are enabled by a uniform gradient of the 

0\ 

external magnetic field. We calculate the probabilities of all unwanted nonresonant 



transitions associated with the flip of each spin with nonresonant frequency and 



£^ . with flips of two spins: one with resonant and one with nonresonant frequencies. 

It is shown that these errors oscillate with changing the gradient of the external 

a: 



magnetic field. Choosing the optimal values of this gradient allows us to decrease 
these errors by 50%. 

PACS numbers: 03.67.Lx, 75.10.Jm 

I. INTRODUCTION 

Even with no interaction with the environment, errors in a quantum computer driven 
by radio-frequency (rf) pulses may be generated by the nonresonant action of these pulses 
on different spins. (See, for example, [l|, Q|.) These errors depend on the value of the de- 
tuning from the resonant (useful) transitions. Unwanted transitions with small detunings 
produce the largest errors. Such transitions, called near-resonant transitions, are associated 
with interactions between the spins in the chain. These interactions are required in order to 
implement the conditional quantum logic. But the same interactions can cause unwanted 
transitions characterized by small detunings from resonance and relatively large probabili- 
ties. If one wants to flip the fcth spin in a particular state and does not want to flip the same 
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spin in the other state with different orientations of the neighboring (k + l)th and (k — l)th 
spins, then first of all one has to suppress these near-resonant transitions. For example, if 
the resonant transition is 

I • • • Ofe-iOfclfc+i ...)—> | . . . Ofc_ilfelfe + i . . .), (1) 

then the (unwanted) near-resonant transitions to be suppressed are 

| . . . Ofc-iOfcOfc+i ...)—*■ | . . . Ofc-ilfcOfc+i . . .), | . . . lfc-iOfclfc+i ...)—> | .. . lfc— ilfclfc+i • • •)• (2) 

A general approach for the complete suppression of near-resonant transitions for an arbitrary 
superposition of quantum sates in a scalable Ising spin quantum computer was developed in 

Q 

Since a wavelength of a rf pulse is much larger than a distance between the spins, the 
pulse affects all spins in the chain. Hence, in order to allow a selective excitation of the 
A;th spin, one must make the Larmor frequency of this particular spin different from the 
Larmor frequencies of all other spins, and choose the frequency of the rf pulse to be equal 
to the transition frequency of this kth spin. In Kane's |2| semiconductor quantum computer 
proposal this can be done by applying a voltage to the kth qubit. The approach used in this 
paper, is to applva permanent magnetic field with a large gradient in the direction of the 
spin chain fl fl Q. Because of the magnetic field gradient, the magnetic field is different 
for different spins, so that each spin has its unique Larmor frequency. 

However, the differences in Larmor frequencies of different spins do not completely solve 
the problem of selective excitations. Even if these differences are large, there are small 
probabilities of excitations of unwanted spins with large detunings. This problem is not 
very important in conventional nuclear magnetic resonance spectroscopy nor for a quantum 
computer with a small number of qubits, since the pulse sequences do not contain a very 
large number of pulses. Unlike this, in a quantum computer with sufficiently large number of 
qubits the implementation of logic operations requires thousands of pulses. The errors caused 
by the nonresonant transitions, being small for a single pulse, in general grow linearly [3] 
with the number of pulses, so that eventually accumulation of this type of error would impose 
severe restrictions on the number of qubits in a nuclear or electron spin quantum computer. 
This problem is especially acute because there are no error correction codes for correction 
of this type of error. 
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In this paper we develop an approach which allows one to decrease significantly (by 
50%) the probabilities of unwanted nonresonant transitions in order to decrease the errors 
caused by the rf pulses of quantum protocols. Our approach to this problem is similar 
to the method used for minimization of errors caused by the near-resonant transitions 
Since the probabilities of unwanted nonresonant transitions oscillate in time with frequencies 
proportional to the differences in the Larmor frequencies between the spins, one can select the 
optimal values of these frequency differences (by applying a specific magnetic field gradient) 
to minimize the nonresonant transitions. In spite of a rather refined character of our model 
(in the sense that this system is complicated for the experimental realization with a solid- 
state device), we believe that our approach provides an important tool for the solution of 
similar problems in many nuclear and electron spin quantum computer proposals. 



II. THE ISING SPIN QUANTUM COMPUTER 

The Hamiltonian for an Ising spin chain placed in an external magnetic field can be 
written in the form, 

H n = -J2 u k P k - 2 J J2 m +1 - % E K exp H K* + ¥>„)] + h.c. } = H + V n (t) . (3) 

fc=0 fc=0 fc=0 

Here h = 1; I k = If ± 1\\ If, If, and 1% are the components of the operator of the kth 
spin 1/2; is the Larmor frequency of the fcth spin; J is the Ising interaction constant; 
fl n is the Rabi frequency (frequency of precession around the resonant transverse field in 
the rotating frame); v n and (p n are, respectively, the frequency and the phase of the nth 
pulse. The Hamiltonian (JHJ) is written for the nth rectangular rf pulse. We assume that the 
Larmor frequency difference 5u> = u>k+i — ojk between neighboring spins is independent of 
the spin's number k. Below we omit the index n which indicates the pulse number. The 
long-range dipole-dipole interaction is suppressed by choosing the angle between the chain 
and the external permanent magnetic field to be equal to the magic angle [jj. 

In the interaction representation, the solution of the Schrodinger equation can be written 
in the form 

*(t) = E c, 9(*)l9> ex pH^*). ( 4 ) 

n 

where E q and \q) are, respectively, the eigenvalues and eigenf unctions of the Hamiltonian 
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H . The expansion coefficients satisfy the system of linear differential equations, 

A(t) = -? E ex p[^(^ - ^ - - (s) 

where a = 1 if E q > E\ and <r = — 1 if E q < E\. The states \l) and \q) in Eq. are related 
by a flip of one fc'th spin, so that the total number of terms in the right-hand side of Eq. 
(|HJ) is equal to the number of qubits, L. 
If the condition 

< J < 5lu (6) 

is satisfied, the pulse effectively affects only one kth spin in the chain Q whose frequency 
uJk is close (near-resonant) or equal (resonant) to the frequency of the pulse, v. In this 
approximation, one has to keep only one term (of the L terms) on the right-hand side of 
Eq. (jHJ and the system of coupled differential equations (JHJ splits into 2 L /2 independent 
pairs of equations 0, 0| of the form 

iC m (t) = -^e-^ pm ^C p (t), (7) 

iLjp — ^ e ^my 1 )- 

Here the states \m) and \p) are related by a flip of the kth spin, A pm = E p — E m — u, and 
we suppose that E p > E m . 

III. NONRESONANT TRANSITIONS 

In this Section, we consider the quantum dynamics which causes the nonresonant effects. 
To simplify the notations, we fix the states \p) and \m), the number of the resonant spin k, 
and define 

5 = E p - E m - v = A pm . (8) 

The solution of Eq. (JJJ) is 



C m {t + t) 



f\r\ 5 /Ar N 
cos T -M-sm - 



e"^ r , (9) 



C p (t + r) = zj sin e iSt °~^. 
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Here to is the time of the beginning of the pulse; A = \J5 2 + fi 2 ; r = t — t is the duration 
of the pulse; and the initial conditions are 

C m (t ) = 1, C p (t ) = 0. (10) 

The solution for the transition from the upper state to the lower state is 

C m (to + r)=i jmn (^-) e - iS ^-^\ 



Xt\ .5 . \t 
cos — - 1- sin — 
\ 2 A I 2 



; 5 



e^ r , (11) 



C p (t + r) 

C m {t ) = 0, C p {t ) = 1. 
Suppose that initially only one state is populated, for example, 

C m (* ) = l, C p (* ) = 0, |m) = |0 4 3 2 1 1 ), |p) = IO4O3I2I1O0) (A: = 2). (12) 

The dynamics of the coefficients C m (t) and C p (t) is defined by Eq. @. Now we calculate the 
probabilities for transitions associated with the flip of a spin with nonresonant frequency, 
and with flips of two spins: one with resonant and one with nonresonant frequencies. We 
consider, for example, the transitions |m) — > |z) and \m) — > \j), where 

\i) = IO4I3O2I1O0), \j) = IO4I3I2I1O0). (13) 

Here k' = 3 is the number of the spin with the nonresonant frequency, and the initial 
conditions are 

a(to) = 0, Cj(t ) = Q. (14) 

The dynamics of these coefficients is defined by Eqs. (jHJ) with two essential terms on the 
right-hand side, 

l c i{t ) = _^ e -*(At-*oq.( f ) _ ^m-E m ^)t-^)\ Cm{t ^ 
l c, J {t) = -^ M -^d{t) - ^e^-^-^-^c^t), (15) 

where 

A = Ej — Ei — is; (16) 

£y > -E^ a = 1 if > E m (and £y > E n ) and cr = -1 if Ei < E m (and Ej < E n ). For 
example, for the states in Eqs. (fT2~j) and (|T5jl one has > E m and a = 1. We characterize 
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the nonresonant transition with (large) detuning D (which in general depends on the kind 
of state and spin number k') defined as 

D = E i -E m -av- 5 -^. (17) 

Then using Eqs. (jHJ), (fTBj) . and (fTTj) one can write 

5 - A 5 - A 
Ei-E m -av = D + — — , E 3 -E v -ov = D — . (18) 

It is convenient to introduce new coefficients 

d(t) = e-^A,, C 3 {t) = jft-^Aj. (19) 

Then one obtains a system of two coupled differential equations for the coefficients Ai(t) 
and Aj(t) with an effective oscillating external force, 

iMt) - jMt) + ^Mt) = ~ei( D ^y-^M Cp (t), (20) 

Ai(to) = 0» M 1 *) = 0, 

where C m {t) and C p {t) are defined by Eq. (J^J). In order to remove the constant phases from 
Eq. (HU we write 

Ai(t) = SP+^-^B^t), Aj(t) = e'lKHB^r). (21) 

The coefficients Bi(r) and Bj(r) satisfy the following system of two coupled differential 
equations: 



iB i (r) + -B i (r) + -B j (r) = f(T) 



iBj { T ) - ^B 3 {r) + % t {r) = g(r), (22) 



B i (0)=0, 5,(0) = 0, 
where the differentiation is performed with respect to time-interval, r, and 



Q 2 



9[T) = "4A 



(23) 
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By differentiation of the system of equations (|22|). one obtains two uncoupled second order 
differential equations, 

Bi(r) + ^U(r) = |/(r) + ^g(r) - if(r), 



and 



Bi(0) = 0, 5,(0) = z| 



(24) 



5,(0) =0, 4(0) = 0, (25) 

where the initial conditions for -Bj(O) and -8,(0) follow from Eq. (J22j) with /(0) = — and 
g(0) = 0. 

The solution for both coefficients Bi(r) and Bj(t) has the form, 



with different coefficients Keeping only first order terms in Q/D, one obtains, 



Rn « J Mr 

±5j(r) = < cos — 

V ; 2D I V 2 



/Ar\ .5 . Xt 
cos — + 1— sm — 
V 2 ) A I 2 



Si(r) 



Q { Q 



-1 < — sm , 

2D I A V 2 



Xt 



JDt 



(26) 



This solution is written for the initial condition ()10|). For the initial conditions given by the 
third equation (fTTj) . the solution is 

B'(t) = -1 -sm — - — sm — e } , 

tK ' 2D 1 A V 2 / A 2 I 



5 



'At* 



cos T -z-sm T 



JDt 



(27) 



Here the coefficients B'^r) and B'-{r) are obtained using Eq. (JTTJl and the properly modified 
relations (|19|) and (|2ip (which result in different common phase factors, unimportant for 
further consideration) . 
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IV. PROBABILITIES OF NONRESONANT TRANSITIONS 

The nonresonant transitions considered in the previous Section generate errors in the 
quantum computer. Our objective is to minimize the probabilities of these transitions (for 



estimations of these probabilities see Refs. p], I&I9J) as much as possible by choosing optimal 
parameters of the model. As follows from Eq. (j2H|), the probability of a nonresonant 
transitions associated with the flip of a spin with nonresonant frequency is proportional 
to (Q/25u) 2 . (See also Refs. j?l If one could suppress these transitions, the total 
probability of error would be associated with flips of two spins with nonresonant frequencies, 
and it would be of order of (Q/25uj) a . Successful solution of this problem would allow one to 
substantially decrease the errors introduced by the rf pulses and to relieve the requirements 
for large magnetic field gradients characterized by 5u. 

In order to estimate the errors given by Eq. (|26jh we will analyze this equation for 
a typical example. Consider the errors generated during implementation of the (useful) 
transition 

\m) -> \p), (28) 

where the states \m) and \p) are defined in Eq. ()12jl . Since the frequency of the external rf 
field is resonant for this transition, then 

v = E v -E m = 35lu, 6 = \ = tt, r = ^ k = 2, (29) 

where k is the number of the spin with the resonant frequency. 
The values 



A k> , \ k , = ^/Al + W (30) 
for nonresonant transitions to the states with one flipped k' spin are 

A fc / = 5, A fc / = A if fcV k ± !• ( 31 ) 

When kl = ±1 the detunings are different, ^ 5, and depend on the orientations of the 

(k + l)th and (k — l)th spins in the state For example, for the states in Eq. (|13|) one 
has Ej — Ei = 35u — 2 J so that 

A 3 = Ej - Ei - v = 5 - 2 J = -2 J. (32) 
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For the detunings associated with the other spins, one has 

A = 0, Ax = 2J A 4 = 0. (33) 



In order to produce resonant transitions (where 5 = 0), the time-interval should be 
t = ir/Q [7r-pulse, see Eq. (J2J)]. In order to suppress possible near- resonant transitions, the 
value of Q should satisfy the 27rA'-condition 



iupj 



V4K 2 - 1 V ! 

For these parameters and for the state \m) in Eq. ()12)1 (for k = 2) the quantities in Eq. (j2 
take the following values: 

5 = 0, X = n, sin i^L j = 1, cos j = 0, 



sin 



sin 



, 1, cos (^-|- ) = 0, for k' ^ k± 1, (35) 
Afc ' T ^ 0, cos^l I . tor = A- ± 1. 



2 / V 2 

Then, after the 7r-pulse the coefficients in Eq. 1)26)1 become 



5f=0, S* =-z^-(l-e-^) for fc^ * ± 1, (36) 



and 



B? = -■£-, fif=i-5_ e ^ for A> = *±1, (37) 

where the upper index of the amplitudes Bf and Bj' indicates that these states are associ- 
ated with the initial state by the flip of the k'th spin. 

The coefficients B? and Bf in Eqs. © and §7\) are calculated for a transition resonant 
for the initial state. If the frequency of the rf pulse is close to resonant (near- resonant), then 

5^0, A = VQ 2 + 4 J 2 , sin (^pj = 0, cos \^-J = 1, 

f)=0, cos(^)=l, for k'^k±l, (38) 



sm 



sin 



1, cos ( -|- ) = 0, for k' = k± 1. 
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FIG. 1: Probability of errors after application of 7r-pulse resonant for the transition (|28|) for states 
and initial conditions denned in Eq. 1)12(1 . The solid line indicates the analytic results calculated 
using Eq. (|39|l .The filled circles are calculated using exact numerical solution 0]. L = 5, 
n/J= 0( A ")/J ss 1/2, K = 2 [see Eq. 

The error for an arbitrary initial state generated by the pulse resonant for the transition 
is the sum of the probabilities of all nonresonant transitions associated with flips of all 
spins with k! ^ k, 



Pk 



n 



n 



2Dk-i 



2D, 



+ 



E 



\2Dy 



(39) 



If k — 1 then the term with D^_i in Eq. (|3*9*j) should be omitted. If k — L — 1, one must 
omit the term with -Dfc+i- In Fig. 1 we plot the total probability of error after application 
of a 7r-pulse resonant for the transition with the states and initial conditions defined in 
Eq. (JT2J) . One can see a good correspondence between the numerical and analytical results. 
The same agreement was observed for other initial states and parameters of the model. The 
parameters used in the simulations can be realized using, for example, phosphorus impurity 
donors in silicon. (For the relation between the numerical and physical parameters, see, for 
example, Ref. 
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FIG. 2: The same as Fig. 1 but for different states. The results are obtained using Eq. (|39|). 



V. MINIMIZATION OF NONRESONANT TRANSITIONS 



From Eq. (J39j) and Fig. 1 one can see that probability of error oscillates as a function of 
the ratio 5uj/VL. As follows from Eq. (p?9*|) the period of oscillations is equal to 2. (In Fig. 
1 the period is twice smaller because of the symmetry of the transition: the probabilities 
oscillate only for the transitions with detunings of the order of 25u.) 

In order to describe these oscillations we continue to analyze the example from the pre- 
vious Section. The values of Dy for the initial state \m) in Eq. (|T2"|) and for k = 2 are 



-25uj -J, Dx = 6oj - J, D 4 = 5u + J, D 5 = 25lu + J. 



(40) 



Under the condition 



|A>, ,D 5 



(41) 



all (two) terms in the sum in Eq. (|39|) vanish, and the probability of error takes its minimum 
value equal to 

pn " - Q 2 - ?■ (42) 

For example, in Fig. 2 for 5uj/VL m q = 194 one has P mm w 2.65 x 10~ 5 . 

We analyzed the nonresonant transitions only for one state. We now discuss the general 
case. In Fig. 2 we plot the probabilities of errors for four randomly selected states generated 
by the pulse, resonant for the transition ((TJ) for L = 20 qubits. For comparison, we also plot 
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the curve from Fig. 1 (solid line) analyzed before. In principle, one could calculate errors 
for an arbitrary number of qubits. However, since the probabilities of flips of distant qubits 
decrease as l/(k — k') 2 , the distant qubits do not contribute significantly to the total error. 

From Fig. 2 one can see that errors for different states take their minimum values at the 
same values of 5u/Q. It is convenient to write the detunings Dkk' (here k is the number 
of spin with the resonant or near-resonant frequency, and k! is the number of spin with 
nonresonant frequency) in the form 

\D k k'\ = \k - k'\Su + a kk >J, (43) 

where a kk i can assume the integer values —2, —1, 0, 1, 2. Choosing the parameters 

^ = 2Q, K = yC, (44) 

where Q ^> 1 and JC are integers, one can write 

Dkk' 



n 



Here we used Eq. The value of Q ^> 1 defines the magnetic field gradient characterized 

by the ratio 5u/Q, and /C > 1 determines the Rabi frequency, or Q/ J at a given value of 
the Ising interaction constant J. 

Under the conditions (jB|) the errors Pkk' associated with flips of k'th spins with k! ^ 
k, k ± 1 (below called the distant spins) in the sum of Eq. (f3T)|) take their minimum values 
equal to 

At the minimums, the transitions associated with the flips of the distant spins are mostly 
suppressed, and the total error ()42)1 is mainly associated with flips of the (k — l)th and the 
(k + l)th spins, 



pmm 



4Q 2 



1 2/C / I \ ( 7T \ 2 _ ( a kk ' 



Q V 32/C 2 / V32/C/ fc /^ fc _ lifc)fe+1 V/c fc' 

(47) 

From this equation one can see that even for minimum value K — 1, the error associated 
with flips of distant spins is at least two orders smaller than the error associated with flips 
of (k + l)th and (k — l)th spins, and this error quickly decreases as K, increases. 
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FIG. 3: The error for different number L of spins in the spin chain, k is chosen in the center of 
the chain when the error takes its maximum value. KL = 1. 



Since the transitions associated with flips of the distant spins are mostly suppressed, for 
the optimal parameters (J44)) the error is practically independent of the number of qubits. 
In Fig. 3 we compare the errors generated by nonresonant transitions for different numbers 
of spins in the spin chain. From Fig. 3 one can see that, in general, the error grows as 
the number of qubits increases, but the error in the minima is almost independent of the 
number of qubits L. 

It is important that the error is associated only with the (k + l)th, fcth, and the (k — l)th 
spins. The probabilities of flips of other spins are small when the condition ()44|) is satisfied. 
In Figs. 4(a,b) the probabilities of error are plotted for different states for the optimal 
value Su/n = 2Q = 200 [Fig. 4(a)], and for another value 5u/Q = 200.314 [Fig. 4(b)]. 
By comparison of Fig. 4(a) and Fig. 4(b), one can see that optimization of the parameters 
results in the relative suppression of all but four transitions, associated with flips of (k — l)th, 
kth, and (k + l)th spins. Since the errors are associated with the definite three spins these 
errors, probably, can be corrected by application of error-correction codes associated with 
these particular spins. 
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FIG. 4: Probabilities of unwanted states created in result of the 7r-pulse resonant to the transition 
Q. Initial state is 1 11010111001) , k = 5, L = 11, K = 5 (O/J » 0.1); (a) = 200 (Q = 100), 

(b) 6oj/tt = 200.314. The results are obtained using exact numerical solution 0]. 

VI. CONCLUSION 



In this paper, we considered the effect of suppression of nonresonant transitions only for 
one pulse. Since a quantum protocol consists of many such pulses, application of our results 
would substantially decrease the rate with which the error accumulates in the register of a 
quantum computer. 

In our model the Larmor frequency Wq + k'Su of the k'th spin increases with increasing 
spin number k' . If k is the number of the spin with resonant frequency, then the moduli of 
detunings 8cu\k — k'\ increase with increasing "distance" \k — k'\, so that the error associated 
with flips of distant spins decreases as l/(k — k') 2 . The situation is different for Kane's 
quantum computer |3|. In this computer the Larmor frequencies are the same for different 
spins, except for one spin whose frequency is changed by application of a voltage to this 
particular spin. Hence, the detunings are the same for all spins in the chain (except for the 
one spin to which the voltage is applied) and equal to 8uj. As a result, the error caused 
by the nonresonant transitions should increase linearly with an increase in the number of 
qubits. Our method would allow one to suppress this linear growth. 
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